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Abstract. In this note we discuss a couple of teciinical issues relevant to solving the Dyson-Schwinger equation for the 
gluon propagator in Landau gauge Yang-Mills theory. In the deep infrared functional methods extract a one-parameter family 
of solutions generically showing a massive behavior referred to as 'decoupling' but also including the so-called 'scaling' 
solution with a conformal infrared behavior as a limiting case. We emphasize that the latter cannot be ruled out by technical 
arguments related to the removal of quadratic divergencies and transversality. 
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Introduction 

The infrared behavior of the Green's functions of Lan- 
dau gauge QCD has been focus of a number of discus- 
sions over the last years. In general, this interest is mo- 
tivated by the desire to understand global properties of 
the theory such as the realization of the global gauge 
charges which is closely related to the description of con- 
finement, as formulated by Kugo and Ojima many years 
ago [7], within the local field theory framework of co- 
variantly gauge-fixed Yang-Mills theory. Clearly, non- 
perturbative approaches such as lattice QCD or func- 
tional methods are necessary to test these ideas. 

From the functional methods, be it Dyson-Schwinger 
equations (DSEs) or the Functional Renormalization 
Group (FRG), it is now known that two types of solu- 
tions exist which differ in their asymptotic infrared be- 
havior One is a 'massive' or 'decoupling' type of solu- 
tion, which is characterized by an infrared finite gluon 
propagator and a ghost propagator with an infrared fi- 
nite dressing function ylS 13.1511^ and the other is the so- 
called 'scaling' solution with unique infrared power laws 
for both propagators S 0> HI 01 and all other Green's 
functions of Landau gauge Yang-Mills theory flo',nT]. 

In terms of the dressing functions G{p^) and Z{p^) of 
the ghost and gluon propagators in Landau gauge 

^cyP) = -p2- ' -^MvlP) = [Otiv - j , 

(1) 

decoupling is characterized by the infrared behavior 

Z{p^) - p^/M^ , G{p^) const. , for 0, (2) 

whereas with scaling in 4 dimensions one has, 

Z(/) G(/)^(/)-^ (3) 



with a positive exponent fc < 1 which under a certain reg- 
ularity assumption on the ghost-gluon vertex [7] results 
to be K- = K-^ = (93 - VTm) /98 w 0.6. 

Both, scaling and the decoupling type of solutions 
together form a one-parameter family, which has been 
obtained analytically and as infrared limits of numerical 
solutions to functional equations. With fixed gluon input 
this family has been found in the ghost DSE in Ref. |3]; 
full numerical solutions for the coupled ghost and gluon 
system of DSEs and FRGEs have been given for the 
whole one-parameter family in Ref. |5]. The decoupling 
results agree quantitatively well with lattice data |12]. 
Scaling has not been observed unambiguously on the 
lattice so far. There is an ongoing effort, however, to 
understand why that is the case and what the differences 
are between the functional continuum methods and those 
commonly used in what is called lattice Landau gauge 

Meanwhile, we would like to reply in this note to a 
claim made by Mike Pennington concerning an apparent 
problem with scaling from continuum arguments alone. 
In his plenary talk at this conference he purported that 
infrared scaling was observed only as a result of an in- 
consistency between the way in which quadratic diver- 
gences are being removed and the transversality of the 
gluon DSE in Landau gauge. In the following we explain 
why this conclusion is itself incorrect. There is no prob- 
lem with quadratic divergences and transversality in the 
scaling results from the functional continuum methods. 



Transversality of the gluon propagator 

In order to understand the matter let us first recall, that 
in covariant gauges there is a Slavnov-Taylor identity for 



the gluon propagator stating that its longitudinal part is 
not modified by interactions, 



d^dyD;Ux-y) = ^5'"'5\x-y) 



ab^Ar 



(4) 



In the Landau gauge limit ^ it then follows 
that the fully dressed gluon propagator remains trans- 
verse which implies that its momentum space Dyson- 
Schwinger equation is of the form, 



(5) 



where Z3 is the gluon renormalization constant, <^^y = 

(5fiv — PfiPv/ p^) the transverse projector, and Tl^v{p'^) 
the self-energy which is then necessarily transverse also. 
In practical calculations, however, the transversality of 
n^v is often difficult to maintain in specific truncations. 
The main sources of transversality violations are thereby 
the following: (i) Numerical solutions of dimension- 
ally regularized integral equations are extremely cum- 
bersome 1 16]. In practice one therefore uses different 
schemes such as momentum subtractions with a hard cut- 
off. While these schemes in general preserve multiplica- 
tive renormalizability of the theory, they violate Eq. (|4|l. 
(ii) Ansaetze for the different vertices in Yl^vip'^), nec- 
essary to close the gluon DSE, may be inadequate to pre- 
serve transversality. In both cases, (i) and (ii), artificial 
longitudinal contributions £^^yL{p^) arise on the right 
hand side of the gluon DSE, 



(6) 



with B^^y — PiiPv/p^- A different but closely related 
problem is the appearance of quadratic divergences in 
n^v Again, these are absent in dimensional regulariza- 
tion, but they occur with the cutoff procedure typically 
used in numerical studies of the gluon DSE i2[l5L ll7ll . 

When discussing the non-perturbative infrared prop- 
erties of the gluon, care must be taken that artifacts like 
non-transversality or the appearance of quadratic diver- 
gencies do not affect the results. This may be particularly 
important when it comes to the discussion of scaling vs 
decoupling, Eqs.©, ([3]l. In his talk, Mike Pennington has 
addressed this problem, claiming that the very appear- 
ance of the scaling solution di]) is such an artifact. His 
argument relies on the assumption that the ghost-gluon 
vertex is essentially bare in the infrared. If one then uses a 
particular method to remove quadratic divergencies (see 
Eq. ([8j below) the scaling solution disappears. 

This problem has been addressed already in Ref. 0]: 
Scaling goes hand-in-hand with the infrared dominance 
of ghosts, i.e. the ghost-loop dominates Y\^y{p^), i.e., 
for p^ 0. Therefore, the ghost-loop must itself be 
transverse in the infrared. In a truncation where the full 



ghost-gluon vertex is replaced by the tree-level vertex 
of standard Faddeev-Popov theory this would require 
ff = 3/4, which is incompatible with the self-consistently 
obtained value fc^. « 0.6 mentioned above. Thus the non- 
perturbative ghost-gluon vertex can not be identical to 
the tree-level one. 

Instead, infrared dominance of ghosts and with that the 
scaling solution essentially requires that the fully dressed 
ghost-gluon vertex becomes itself transverse with re- 
spect to the gluon momentum in the infrared [7]. Under 
the mild additional regularity assumption one then ob- 
tains the self-consistent value k~Kc~ 0.6. Other values 
1 /2 < K < 1 are possible when this condition is relaxed 
1 7]. Note, that all these analytic results were obtained in 
dimensional regularisation. They are thus unaffected by 
any technical difficulties that might arise with quadratic 
divergences in numerical investigations. 

One method to remove quadratic divergences in the 
gluon DSE goes back to Brown and Pennington ifisll . 
They decomposed the gluon self-energy into 



n^vb^) = d^vF\{p^)- p^iPvFiip^) 



(7) 



and noted that quadratic divergences can only occur in 
the term proportional to 8fiv Consequently, they sug- 
gested to project out this contribution. This is don e by 
contracting the gluon-DSE with a general projector ilTIl 



c 



PjiPv 



(8) 



and setting C, — A in d = A dimensions. This removes 
quadratic divergencies and leads to a dependence of 
Z[p'^) on F2{p^) alone. Alas, let us see what one ob- 
tains for general ^ upon contracting the gluon-DSE in 
the form of Eqs. Q and ^ with iP^y- 



Z{p^) 



= p2z3+n(p2 



Li.p'-) 



(9) 



Clearly, if the gluon DSE is transverse, Llp^) — 0, and 
the resulting Eq. (|9]l is independent of the parameter 1^ 
and therefore also free from quadratic divergences. Con- 
versely, the required (^-independence provides a valuable 
test of specific truncations in numerical studies. 

Two transverse truncations passing this test have 
been constructed explicitly in the literature, and both 
of them rely on a nontrivial ghost-gluon vertex thus 
dismissing Mike's criticism of the inconsistency that 
arises with the bare one. One is based on the Pinch- 
Technique/Background Field Method, see fl^, for 
which the truncation of Ref |2] has also been analyzed 
analytically in dimensional regularizationQ The other 



' Note, however, that the numerical solutions presented in Ref. @] have 
been obtained using a hard UV-cutoff and therefore need to be checked 
separately for transversality. 



approach is the one used in Ref. lHj]. There, Ansaetze for 
the ghost-gluon and three-gluon vertices have been con- 
structed expHcitly such that the gluon-DSE is free from 
quadratic divergencies and the residual contributions 
to L{p^) are minimized. The ghost-gluon vertex used 
there reduces in the infrared to the minimally dressed 
transverse one proposed in It contains a bare part 
and a nontrivial longitudinal part that serves to cancel 
all contributions to L{p^), thus explicitly establishing 
exact transversality of the ghost-loop in the infrared. 
Numerically, the remaining contributions to L{p^) for 
all momenta stay well below one percent of n(p^). 
In particular L{p^) — > for p-^ 0. This therefore es- 
tablishes unambiguously that the DSE results of fsl] are 
not affected by transversality violating artifacts: there 
is a one-parameter family of solutions in the infrared 
including decoupling and scaling solutions. 

The earlier truncation scheme from Ref. |17], which 
was particularly criticized by Mike Pennington in his ple- 
nary talk, was indeed less far developed and not mani- 
festly transverse. The artificial longitudinal term L{p~) 
in Eq. (|9]l was removed by transverse projection, using 
C — 1 rather than ^ = 4 as suggested by Brown and Pen- 
nington before, which meant that quadratic divergences 
had to be dealt with in a different manner Even though 
this procedure might seem less elegant, there were good 
reasons at the time for that given in iQ] . Before we briefly 
explain those, it can now be verified a posteriori, by com- 
paring the earlier results from [17] to those of Ref. jstl, 
to verify that they are almost identical (for ^ = 1). 

The reason why transverse projection onto Tl{p^) in 
Eq. ^ must be done |7] rather than ^ ~4 projection onto 
F2{q^) in (|7]i, in a not manifestly transverse truncation 
where this can make a difference, relates to an ambigu- 
ity in the Landau gauge: the existence of a second gauge 
parameter which interpolates between standard Faddeev- 
Popov theory and the ghost/anti-ghost symmetric Curci- 
Ferrari gauges. In Landau gauge there is no distinction, 
however, and any quantity depending on this parameter is 
thus ambiguous. In [TJ it was shown that /^2(f ^) is such 
a quantity while n(/7^) is not. So we have understood 
over the years that it is the original proposal from fis!] 
to avoid quadratic divergences which is ambiguous, un- 
fortunately. Luckily, however, the problem is now com- 
pletely solved with the manifestly transverse truncations 
that were developed since then. 

Finally, we emphasize that the problem with quadratic 
divergences is completely unheard of in the FRGEs 
which are finite by construction. Yet, one still finds the 
full one-parameter family of scaling and decoupling so- 
lutions [9, 5]. Apart from the extreme infrared, where this 
parameter matters, these solutions afl agree and are fur- 
thermore in almost perfect agreement with lattice data. 



Conclusions 

In this note we have pointed out that a claim made by 
Mike Pennington in his plenary talk is not correct. He 
purported that the infrared scaling solution for the ghost 
and gluon propagators of Landau gauge Yang-Mills the- 
ory is a mere artifact of a technical inaccuracy in the 
treatment of the gluon DSE. This is based an the overly 
simplistic truncation which cannot be combined with 
an earlier proposal to remove artificial quadratic diver- 
gences in an untlraviolet cutoff regularization. This is- 
sue has been completely solved over the years, however. 
The infrared-scaling solution for the Yang-Mills sector 
of QCD does not have fundamental problems related to 
quaderatic divergences or transversality. 
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